{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}C:\Users\jcjo228\Dropbox\DeterResponse\LeedsJohnsonJOPrep\LeedsJohnsonJOPrep.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}16 Aug 2015, 15:52:01

{com}. do "C:\Users\jcjo228\Dropbox\DeterResponse\LeedsJohnsonJOPrep\LeedsJohnsonJOPrep.do"
{txt}
{com}. * This file runs the analysis discussed in "Theory, Data, and Deterrence:  A Response to Kenwick, Vasquez, and Powers"
. * Brett Ashley Leeds and Jesse C. Johnson
. * 8/16/15
. 
. **********************************
. *** Code for Empirical Results ***
. **********************************
. 
. set more off
{txt}
{com}. cd "C:\Users\jcjo228\Dropbox\DeterResponse\LeedsJohnsonJOPrep"
{res}C:\Users\jcjo228\Dropbox\DeterResponse\LeedsJohnsonJOPrep
{txt}
{com}. use "LeedsJohnsonJOPrep.dta", clear
{txt}
{com}. 
. * We replicate the JL FPA results (Table A1 and Figure A1)
. 
.         estsimp probit dispute ptargdef pchalally pchaloff pchalneu ln_distance capprop jdem s_un_glo peaceyrs peaceyrs2 peaceyrs3 

{txt}Iteration 0:   log likelihood = {res}-16773.772
{txt}Iteration 1:   log likelihood = {res}-14885.021
{txt}Iteration 2:   log likelihood = {res}-14127.643
{txt}Iteration 3:   log likelihood = {res}-14089.817
{txt}Iteration 4:   log likelihood = {res}-14089.525
{txt}Iteration 5:   log likelihood = {res}-14089.525

{txt}Probit regression                                 Number of obs   = {res}   1077992
                                                  {txt}LR chi2({res}11{txt})     = {res}   5368.49
                                                  {txt}Prob > chi2     = {res}    0.0000
{txt}Log likelihood = {res}-14089.525                       {txt}Pseudo R2       = {res}    0.1600

{txt}{hline 13}{c TT}{hline 64}
     dispute {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
    ptargdef {c |}  {res}-.0615538   .0154694    -3.98   0.000    -.0918733   -.0312342
   {txt}pchalally {c |}  {res} -.067629   .0156073    -4.33   0.000    -.0982187   -.0370394
    {txt}pchaloff {c |}  {res} .2762056   .0286835     9.63   0.000      .219987    .3324243
    {txt}pchalneu {c |}  {res} .3146708   .0251969    12.49   0.000     .2652859    .3640557
 {txt}ln_distance {c |}  {res}-.3836792   .0064471   -59.51   0.000    -.3963154   -.3710431
     {txt}capprop {c |}  {res} .1611667   .0224281     7.19   0.000     .1172084    .2051249
        {txt}jdem {c |}  {res}-.1545372   .0266195    -5.81   0.000    -.2067105   -.1023638
    {txt}s_un_glo {c |}  {res}-.4553676   .0412643   -11.04   0.000    -.5362441   -.3744911
    {txt}peaceyrs {c |}  {res}-.0420438   .0013683   -30.73   0.000    -.0447257   -.0393619
   {txt}peaceyrs2 {c |}  {res} .0006635    .000028    23.66   0.000     .0006085    .0007184
   {txt}peaceyrs3 {c |}  {res}-2.78e-06   1.48e-07   -18.82   0.000    -3.07e-06   -2.49e-06
       {txt}_cons {c |}  {res} .8266646   .0675454    12.24   0.000     .6942781    .9590511
{txt}{hline 13}{c BT}{hline 64}

{res}Simulating main parameters.  Please wait....
% of simulations completed: 8% 16% 25% 33% 41% 50% 58% 66% 75% 83% 91% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12
{txt}
{com}. 
.         setx mean if ptargdef==0
{txt}
{com}.         setx ptargdef 0
{txt}
{com}. 
.         simqi, prval(1) genpr (nodef)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0009437      .000041     .0008657    .0010277

Simqi generated the following new variable(s): nodef
{txt}
{com}.         
.         _pctile nodef, p(2.5,97.5)
{txt}
{com}.         replace nodef =. if nodef < r(r1) | nodef > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize nodef

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}nodef {c |}{res}        950    .0009436     .000036   .0008657   .0010275
{txt}
{com}. 
.         setx mean if ptargdef==1
{txt}
{com}.         setx ptargdef 1
{txt}
{com}. 
.         simqi, prval(1) genpr (def)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0007111     .0000311     .0006496    .0007707

Simqi generated the following new variable(s): def
{txt}
{com}.         
.         _pctile def, p(2.5,97.5)
{txt}
{com}.         replace def =. if def < r(r1) | def > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize def

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 9}def {c |}{res}        950    .0007109     .000027   .0006502   .0007703
{txt}
{com}. 
.         gen def_change = ((def-nodef)/nodef)*100
{txt}(1,077,084 missing values generated)

{com}.         sum def_change

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 2}def_change {c |}{res}        908   -24.56266    3.412665  -34.95952  -14.81386
{txt}
{com}.         
.         twoway hist def, frequency color(gray 8) ||hist nodef, frequency color(black) legend(order(1 "Target with Defense Pact" 2 "Target without Defense Pact")) ysc(on) xlabel(0 .0005 .001) xtitle("Probability of Dispute Initiation") ylabel(0 20 40 60) ytitle(Frequency) scheme(s1mono) 
{res}{txt}
{com}.         capture drop b1-b12 nodef def def_change
{txt}
{com}. 
. * We replicate the Figure 1A KVP model (Table A2 Column 1)
. 
. use "KVPJOPrep.dta", clear
{txt}
{com}.         
.         cem jdem(#0) pchaloff(#0) pchalneu(#0) contig(#0) ccode1_maj(#0) ccode2_maj(#0) rivalry_bef(#0) disputes_before (#0) ptargoff(#0) ptargneu(#0), tr(longitudinal_treatment)
{res}
{txt}Matching Summary:
-----------------
Number of strata: {res}311
{txt}Number of matched strata: {res}89

           {txt}    0      1
      All  {res}38732   6989
{txt}  Matched  {res}38059   6870
{txt}Unmatched  {res}  673    119


{txt}Multivariate L1 distance: {res}2.238e-13

{txt}Univariate imbalance:

                       L1      mean       min       25%       50%       75%       max
           jdem  {res} 5.3e-14  -1.4e-14         0         0         0         0         0
{txt}       pchaloff  {res} 5.5e-14  -1.0e-14         0         0         0         0         0
{txt}       pchalneu  {res} 8.5e-14  -1.4e-14         0         0         0         0         0
{txt}         contig  {res} 4.1e-14  -8.4e-15         0         0         0         0         0
{txt}     ccode1_maj  {res} 8.1e-14  -1.6e-14         0         0         0         0         0
{txt}     ccode2_maj  {res} 1.1e-13  -1.8e-14         0         0         0         0         0
{txt}    rivalry_bef  {res} 1.1e-14  -1.4e-14         0         0         0         0         0
{txt}disputes_before  {res} 3.9e-15  -1.7e-15         0         0         0         0         0
{txt}       ptargoff  {res} 4.7e-14  -9.2e-15         0         0         0         0         0
{txt}       ptargneu  {res} 6.2e-14  -9.4e-15         0         0         0         0         0
{txt}
{com}. 
.         logit anydispaft longitudinal_treatment [iweight=cem_weights]

{res}{txt}Iteration 0:{space 3}log likelihood = {res:-2250.7794}  
Iteration 1:{space 3}log likelihood = {res:-2247.1993}  
Iteration 2:{space 3}log likelihood = {res:-2247.1152}  
Iteration 3:{space 3}log likelihood = {res:-2247.1151}  
{res}
{txt}Logistic regression{col 49}Number of obs{col 67}= {res}    44,929
{txt}{col 49}LR chi2({res}1{txt}){col 67}= {res}      7.33
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0068
{txt}Log likelihood = {res}-2247.1151{txt}{col 49}Pseudo R2{col 67}= {res}    0.0016

{txt}{hline 23}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}            anydispaft{col 24}{c |}      Coef.{col 36}   Std. Err.{col 48}      z{col 56}   P>|z|{col 64}     [95% Con{col 77}f. Interval]
{hline 23}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
longitudinal_treatment {c |}{col 24}{res}{space 2} .3532531{col 36}{space 2}  .125996{col 47}{space 1}    2.80{col 56}{space 3}0.005{col 64}{space 4} .1063054{col 77}{space 3} .6002007
{txt}{space 17}_cons {c |}{col 24}{res}{space 2}-4.794432{col 36}{space 2} .0568131{col 47}{space 1}  -84.39{col 56}{space 3}0.000{col 64}{space 4}-4.905784{col 77}{space 3}-4.683081
{txt}{hline 23}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. 
. * We estimate the KVP model without matching (Table A2 Column 2)
. 
.         logit anydispaft longitudinal_treatment pchaloff pchalneu ptargoff ptargneu contig ccode1_maj ccode2_maj jdem rivalry_bef disputes_before

{res}{txt}Iteration 0:{space 3}log likelihood = {res:-2641.8481}  
Iteration 1:{space 3}log likelihood = {res: -2389.726}  
Iteration 2:{space 3}log likelihood = {res:-2109.8963}  
Iteration 3:{space 3}log likelihood = {res:-1966.2823}  
Iteration 4:{space 3}log likelihood = {res:-1941.3921}  
Iteration 5:{space 3}log likelihood = {res:-1941.1694}  
Iteration 6:{space 3}log likelihood = {res:-1941.1692}  
Iteration 7:{space 3}log likelihood = {res:-1941.1692}  
{res}
{txt}Logistic regression{col 49}Number of obs{col 67}= {res}    45,721
{txt}{col 49}LR chi2({res}11{txt}){col 67}= {res}   1401.36
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log likelihood = {res}-1941.1692{txt}{col 49}Pseudo R2{col 67}= {res}    0.2652

{txt}{hline 23}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}            anydispaft{col 24}{c |}      Coef.{col 36}   Std. Err.{col 48}      z{col 56}   P>|z|{col 64}     [95% Con{col 77}f. Interval]
{hline 23}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
longitudinal_treatment {c |}{col 24}{res}{space 2} .3448984{col 36}{space 2} .1320979{col 47}{space 1}    2.61{col 56}{space 3}0.009{col 64}{space 4} .0859913{col 77}{space 3} .6038056
{txt}{space 14}pchaloff {c |}{col 24}{res}{space 2}-.2186613{col 36}{space 2} .2093337{col 47}{space 1}   -1.04{col 56}{space 3}0.296{col 64}{space 4}-.6289479{col 77}{space 3} .1916252
{txt}{space 14}pchalneu {c |}{col 24}{res}{space 2}-.1660439{col 36}{space 2} .1873464{col 47}{space 1}   -0.89{col 56}{space 3}0.375{col 64}{space 4} -.533236{col 77}{space 3} .2011482
{txt}{space 14}ptargoff {c |}{col 24}{res}{space 2} .3450545{col 36}{space 2} .2826623{col 47}{space 1}    1.22{col 56}{space 3}0.222{col 64}{space 4}-.2089534{col 77}{space 3} .8990625
{txt}{space 14}ptargneu {c |}{col 24}{res}{space 2}-.2284014{col 36}{space 2} .2175618{col 47}{space 1}   -1.05{col 56}{space 3}0.294{col 64}{space 4}-.6548146{col 77}{space 3} .1980118
{txt}{space 16}contig {c |}{col 24}{res}{space 2} 2.151516{col 36}{space 2} .1208096{col 47}{space 1}   17.81{col 56}{space 3}0.000{col 64}{space 4} 1.914733{col 77}{space 3} 2.388298
{txt}{space 12}ccode1_maj {c |}{col 24}{res}{space 2}  1.74637{col 36}{space 2} .1136849{col 47}{space 1}   15.36{col 56}{space 3}0.000{col 64}{space 4} 1.523551{col 77}{space 3} 1.969188
{txt}{space 12}ccode2_maj {c |}{col 24}{res}{space 2} .6163965{col 36}{space 2} .1549838{col 47}{space 1}    3.98{col 56}{space 3}0.000{col 64}{space 4} .3126339{col 77}{space 3} .9201592
{txt}{space 18}jdem {c |}{col 24}{res}{space 2}-1.101932{col 36}{space 2} .2864315{col 47}{space 1}   -3.85{col 56}{space 3}0.000{col 64}{space 4}-1.663327{col 77}{space 3}-.5405365
{txt}{space 11}rivalry_bef {c |}{col 24}{res}{space 2} .2349277{col 36}{space 2} .0300508{col 47}{space 1}    7.82{col 56}{space 3}0.000{col 64}{space 4} .1760292{col 77}{space 3} .2938262
{txt}{space 7}disputes_before {c |}{col 24}{res}{space 2} 1.259876{col 36}{space 2} .0987341{col 47}{space 1}   12.76{col 56}{space 3}0.000{col 64}{space 4} 1.066361{col 77}{space 3} 1.453391
{txt}{space 17}_cons {c |}{col 24}{res}{space 2}-5.566784{col 36}{space 2} .0809182{col 47}{space 1}  -68.80{col 56}{space 3}0.000{col 64}{space 4}-5.725381{col 77}{space 3}-5.408187
{txt}{hline 23}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}.                 
. * We estimate the JL model with several different mataching strategies (Tables A3-6 & Figure A2)
. 
. * Code for the matching procedure is included below 
. 
. use "LeedsJohnsonJOPrep.dta", clear
{txt}
{com}. 
. * 1
.         
.         probit dispute ptargdef [iweight=cem_weights_1]

{res}{txt}Iteration 0:{space 3}log likelihood = {res:-22699.688}  
Iteration 1:{space 3}log likelihood = {res:-22274.622}  
Iteration 2:{space 3}log likelihood = {res:-22260.374}  
Iteration 3:{space 3}log likelihood = {res:-22260.354}  
Iteration 4:{space 3}log likelihood = {res:-22260.354}  
{res}
{txt}Probit regression{col 49}Number of obs{col 67}= {res} 1,029,856
{txt}{col 49}LR chi2({res}1{txt}){col 67}= {res}    878.67
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log likelihood = {res}-22260.354{txt}{col 49}Pseudo R2{col 67}= {res}    0.0194

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}     dispute{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}ptargdef {c |}{col 14}{res}{space 2}-.3500406{col 26}{space 2} .0122705{col 37}{space 1}  -28.53{col 46}{space 3}0.000{col 54}{space 4}-.3740902{col 67}{space 3}-.3259909
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-2.572365{col 26}{space 2} .0070133{col 37}{space 1} -366.78{col 46}{space 3}0.000{col 54}{space 4}-2.586111{col 67}{space 3} -2.55862
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}.         estsimp probit dispute ptargdef pchalally pchaloff pchalneu ln_distance capprop jdem s_un_glo peaceyrs peaceyrs2 peaceyrs3 [iweight=cem_weights_1]

{txt}Iteration 0:   log likelihood = {res}-22699.688
{txt}Iteration 1:   log likelihood = {res}-20291.668
{txt}Iteration 2:   log likelihood = {res}-19716.583
{txt}Iteration 3:   log likelihood = {res}-19694.994
{txt}Iteration 4:   log likelihood = {res}-19694.884
{txt}Iteration 5:   log likelihood = {res}-19694.884

{txt}Probit regression                                 Number of obs   = {res}   1029856
                                                  {txt}LR chi2({res}11{txt})     = {res}   6009.61
                                                  {txt}Prob > chi2     = {res}    0.0000
{txt}Log likelihood = {res}-19694.884                       {txt}Pseudo R2       = {res}    0.1324

{txt}{hline 13}{c TT}{hline 64}
     dispute {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
    ptargdef {c |}  {res}-.3778247   .0133992   -28.20   0.000    -.4040867   -.3515628
   {txt}pchalally {c |}  {res}-.2908032    .013468   -21.59   0.000    -.3172001   -.2644063
    {txt}pchaloff {c |}  {res} .1132444   .0342688     3.30   0.001     .0460788    .1804101
    {txt}pchalneu {c |}  {res} .4350363   .0210349    20.68   0.000     .3938088    .4762639
 {txt}ln_distance {c |}  {res}-.3525037   .0062435   -56.46   0.000    -.3647408   -.3402665
     {txt}capprop {c |}  {res} .2008403   .0185935    10.80   0.000     .1643977    .2372829
        {txt}jdem {c |}  {res}-.2900719   .0290647    -9.98   0.000    -.3470376   -.2331062
    {txt}s_un_glo {c |}  {res}-.8610919   .0375298   -22.94   0.000     -.934649   -.7875348
    {txt}peaceyrs {c |}  {res}-.0347873    .001123   -30.98   0.000    -.0369883   -.0325862
   {txt}peaceyrs2 {c |}  {res} .0005587   .0000219    25.55   0.000     .0005159    .0006016
   {txt}peaceyrs3 {c |}  {res}-2.27e-06   1.12e-07   -20.25   0.000    -2.49e-06   -2.05e-06
       {txt}_cons {c |}  {res} 1.235967    .067546    18.30   0.000     1.103579    1.368354
{txt}{hline 13}{c BT}{hline 64}

{res}Simulating main parameters.  Please wait....
% of simulations completed: 8% 16% 25% 33% 41% 50% 58% 66% 75% 83% 91% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12
{txt}
{com}.         
.         setx mean if ptargdef==0
{txt}
{com}.         setx ptargdef 0
{txt}
{com}. 
.         simqi, prval(1) genpr (nodef)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0025968      .000078     .0024449    .0027548

Simqi generated the following new variable(s): nodef
{txt}
{com}.         
.         _pctile nodef, p(2.5,97.5)
{txt}
{com}.         replace nodef =. if nodef < r(r1) | nodef > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize nodef

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}nodef {c |}{res}        950    .0025965    .0000678    .002448   .0027543
{txt}
{com}. 
.         setx mean if ptargdef==1
{txt}
{com}.         setx ptargdef 1
{txt}
{com}. 
.         simqi, prval(1) genpr (def)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0008164     .0000342      .000749    .0008833

Simqi generated the following new variable(s): def
{txt}
{com}.         
.         _pctile def, p(2.5,97.5)
{txt}
{com}.         replace def =. if def < r(r1) | def > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize def

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 9}def {c |}{res}        950    .0008163    .0000301   .0007491   .0008824
{txt}
{com}. 
.         gen def_change = ((def-nodef)/nodef)*100
{txt}(1,077,089 missing values generated)

{com}.         sum def_change

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 2}def_change {c |}{res}        903   -68.55846    1.240638  -71.86987  -65.17054
{txt}
{com}.         
.         twoway hist def, frequency color(gray 8) ||hist nodef, frequency color(black) legend(order(1 "Target with Defense Pact" 2 "Target without Defense Pact")) ysc(on) xlabel(0 .0005 .001 .0015 .002 .0025 .003) xtitle("Probability of Dispute Initiation") ylabel(0 20 40 60) ytitle(Frequency) scheme(s1mono) 
{res}{txt}
{com}.         capture drop b1-b12 nodef def def_change
{txt}
{com}.         
. * 2
.                 
.         probit dispute ptargdef [iweight=cem_weights_2]

{res}{txt}Iteration 0:{space 3}log likelihood = {res:-14188.034}  
Iteration 1:{space 3}log likelihood = {res:-14009.642}  
Iteration 2:{space 3}log likelihood = {res:-14005.184}  
Iteration 3:{space 3}log likelihood = {res:-14005.179}  
Iteration 4:{space 3}log likelihood = {res:-14005.179}  
{res}
{txt}Probit regression{col 49}Number of obs{col 67}= {res}   776,717
{txt}{col 49}LR chi2({res}1{txt}){col 67}= {res}    365.71
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log likelihood = {res}-14005.179{txt}{col 49}Pseudo R2{col 67}= {res}    0.0129

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}     dispute{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}ptargdef {c |}{col 14}{res}{space 2}-.2876126{col 26}{space 2} .0156172{col 37}{space 1}  -18.42{col 46}{space 3}0.000{col 54}{space 4}-.3182218{col 67}{space 3}-.2570034
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-2.680231{col 26}{space 2} .0086674{col 37}{space 1} -309.23{col 46}{space 3}0.000{col 54}{space 4}-2.697219{col 67}{space 3}-2.663243
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}.         estsimp probit dispute ptargdef pchalally pchaloff pchalneu ln_distance capprop jdem s_un_glo peaceyrs peaceyrs2 peaceyrs3 [iweight=cem_weights_2]

{txt}Iteration 0:   log likelihood = {res}-14188.034
{txt}Iteration 1:   log likelihood = {res}-12885.778
{txt}Iteration 2:   log likelihood = {res}-12584.738
{txt}Iteration 3:   log likelihood = {res}-12573.372
{txt}Iteration 4:   log likelihood = {res}-12573.267
{txt}Iteration 5:   log likelihood = {res}-12573.267

{txt}Probit regression                                 Number of obs   = {res}    776717
                                                  {txt}LR chi2({res}11{txt})     = {res}   3229.53
                                                  {txt}Prob > chi2     = {res}    0.0000
{txt}Log likelihood = {res}-12573.267                       {txt}Pseudo R2       = {res}    0.1138

{txt}{hline 13}{c TT}{hline 64}
     dispute {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
    ptargdef {c |}  {res}-.3270528   .0169475   -19.30   0.000    -.3602693   -.2938363
   {txt}pchalally {c |}  {res}-.3168143   .0174572   -18.15   0.000    -.3510297   -.2825988
    {txt}pchaloff {c |}  {res} .2285158      .0515     4.44   0.000     .1275777    .3294539
    {txt}pchalneu {c |}  {res} .2728394    .040768     6.69   0.000     .1929356    .3527431
 {txt}ln_distance {c |}  {res}-.4069725   .0092704   -43.90   0.000    -.4251421   -.3888028
     {txt}capprop {c |}  {res} .1082951   .0226967     4.77   0.000     .0638104    .1527798
        {txt}jdem {c |}  {res}-.0735119    .038555    -1.91   0.057    -.1490783    .0020546
    {txt}s_un_glo {c |}  {res}-1.273791   .0565725   -22.52   0.000    -1.384671   -1.162911
    {txt}peaceyrs {c |}  {res}-.0385582   .0019015   -20.28   0.000    -.0422852   -.0348313
   {txt}peaceyrs2 {c |}  {res} .0008222    .000055    14.96   0.000     .0007145    .0009299
   {txt}peaceyrs3 {c |}  {res}-4.78e-06   4.06e-07   -11.76   0.000    -5.57e-06   -3.98e-06
       {txt}_cons {c |}  {res} 1.970937   .1052858    18.72   0.000      1.76458    2.177293
{txt}{hline 13}{c BT}{hline 64}
Note: 154 failures and 0 successes completely determined.

{res}Simulating main parameters.  Please wait....
% of simulations completed: 8% 16% 25% 33% 41% 50% 58% 66% 75% 83% 91% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12
{txt}
{com}. 
.         setx mean if ptargdef==0
{txt}
{com}.         setx ptargdef 0
{txt}
{com}. 
.         simqi, prval(1) genpr (nodef)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0021325     .0000765     .0019929    .0022932

Simqi generated the following new variable(s): nodef
{txt}
{com}.         
.         _pctile nodef, p(2.5,97.5)
{txt}
{com}.         replace nodef =. if nodef < r(r1) | nodef > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize nodef

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}nodef {c |}{res}        950    .0021316     .000066    .001993   .0022929
{txt}
{com}. 
.         setx mean if ptargdef==1
{txt}
{com}.         setx ptargdef 1
{txt}
{com}. 
.         simqi, prval(1) genpr (def)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res}  .000727     .0000412     .0006491    .0008081

Simqi generated the following new variable(s): def
{txt}
{com}.         
.         _pctile def, p(2.5,97.5)
{txt}
{com}.         replace def =. if def < r(r1) | def > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize def

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 9}def {c |}{res}        950    .0007269     .000036     .00065   .0008075
{txt}
{com}. 
.         gen def_change = ((def-nodef)/nodef)*100
{txt}(1,077,086 missing values generated)

{com}.         sum def_change

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 2}def_change {c |}{res}        906   -65.89728    1.774752  -70.37693   -61.3022
{txt}
{com}.         
.         twoway hist def, frequency color(gray 8) ||hist nodef, frequency color(black) legend(off) ysc(on) xlabel(0 .0005 .001 .0015 .002 .0025 .003) xtitle("Probability of Dispute Initiation") ylabel(0 20 40 60) ytitle(Frequency) scheme(s1mono) 
{res}{txt}
{com}.         capture drop b1-b12 nodef def def_change
{txt}
{com}.         
. * 3
. 
.         probit dispute ptargdef [iweight=cem_weights_3]

{res}{txt}Iteration 0:{space 3}log likelihood = {res:-16679.912}  
Iteration 1:{space 3}log likelihood = {res: -16458.14}  
Iteration 2:{space 3}log likelihood = {res: -16452.48}  
Iteration 3:{space 3}log likelihood = {res:-16452.474}  
Iteration 4:{space 3}log likelihood = {res:-16452.474}  
{res}
{txt}Probit regression{col 49}Number of obs{col 67}= {res}   866,911
{txt}{col 49}LR chi2({res}1{txt}){col 67}= {res}    454.88
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log likelihood = {res}-16452.474{txt}{col 49}Pseudo R2{col 67}= {res}    0.0136

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}     dispute{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}ptargdef {c |}{col 14}{res}{space 2}-.2947934{col 26}{space 2} .0143338{col 37}{space 1}  -20.57{col 46}{space 3}0.000{col 54}{space 4}-.3228872{col 67}{space 3}-.2666997
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-2.654946{col 26}{space 2} .0080542{col 37}{space 1} -329.64{col 46}{space 3}0.000{col 54}{space 4}-2.670732{col 67}{space 3} -2.63916
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}.         estsimp probit dispute ptargdef pchalally pchaloff pchalneu ln_distance capprop jdem s_un_glo peaceyrs peaceyrs2 peaceyrs3 [iweight=cem_weights_3]

{txt}Iteration 0:   log likelihood = {res}-16679.912
{txt}Iteration 1:   log likelihood = {res}-15018.123
{txt}Iteration 2:   log likelihood = {res}-14620.022
{txt}Iteration 3:   log likelihood = {res}-14605.702
{txt}Iteration 4:   log likelihood = {res}-14605.644
{txt}Iteration 5:   log likelihood = {res}-14605.644

{txt}Probit regression                                 Number of obs   = {res}    866911
                                                  {txt}LR chi2({res}11{txt})     = {res}   4148.54
                                                  {txt}Prob > chi2     = {res}    0.0000
{txt}Log likelihood = {res}-14605.644                       {txt}Pseudo R2       = {res}    0.1244

{txt}{hline 13}{c TT}{hline 64}
     dispute {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
    ptargdef {c |}  {res}-.3290046    .015613   -21.07   0.000    -.3596055   -.2984037
   {txt}pchalally {c |}  {res} -.370414   .0161382   -22.95   0.000    -.4020442   -.3387838
    {txt}pchaloff {c |}  {res} .3363215   .0389109     8.64   0.000     .2600576    .4125855
    {txt}pchalneu {c |}  {res} .4067491   .0296501    13.72   0.000      .348636    .4648623
 {txt}ln_distance {c |}  {res}-.3918727   .0081809   -47.90   0.000     -.407907   -.3758384
     {txt}capprop {c |}  {res} .2065965   .0213093     9.70   0.000     .1648312    .2483619
        {txt}jdem {c |}  {res}-.2284669   .0396517    -5.76   0.000    -.3061828   -.1507511
    {txt}s_un_glo {c |}  {res} -1.09811   .0497421   -22.08   0.000    -1.195603   -1.000617
    {txt}peaceyrs {c |}  {res}-.0306984   .0015516   -19.79   0.000    -.0337394   -.0276574
   {txt}peaceyrs2 {c |}  {res} .0005131   .0000386    13.29   0.000     .0004375    .0005888
   {txt}peaceyrs3 {c |}  {res}-2.32e-06   2.42e-07    -9.57   0.000    -2.79e-06   -1.84e-06
       {txt}_cons {c |}  {res} 1.669443    .091223    18.30   0.000     1.490649    1.848236
{txt}{hline 13}{c BT}{hline 64}

{res}Simulating main parameters.  Please wait....
% of simulations completed: 8% 16% 25% 33% 41% 50% 58% 66% 75% 83% 91% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12
{txt}
{com}. 
.         setx mean if ptargdef==0
{txt}
{com}.         setx ptargdef 0
{txt}
{com}. 
.         simqi, prval(1) genpr (nodef)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0021514     .0000726     .0020093    .0022904

Simqi generated the following new variable(s): nodef
{txt}
{com}.         
.         _pctile nodef, p(2.5,97.5)
{txt}
{com}.         replace nodef =. if nodef < r(r1) | nodef > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize nodef

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}nodef {c |}{res}        950    .0021509    .0000635   .0020104     .00229
{txt}
{com}. 
.         setx mean if ptargdef==1
{txt}
{com}.         setx ptargdef 1
{txt}
{com}. 
.         simqi, prval(1) genpr (def)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0007291     .0000374     .0006547     .000806

Simqi generated the following new variable(s): def
{txt}
{com}.         
.         _pctile def, p(2.5,97.5)
{txt}
{com}.         replace def =. if def < r(r1) | def > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize def

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 9}def {c |}{res}        950     .000729    .0000326   .0006549   .0008059
{txt}
{com}. 
.         gen def_change = ((def-nodef)/nodef)*100
{txt}(1,077,088 missing values generated)

{com}.         sum def_change

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 2}def_change {c |}{res}        904   -66.08475    1.676513  -70.73731  -60.68737
{txt}
{com}.         
.         twoway hist def, frequency color(gray 8) ||hist nodef, frequency color(black) legend(off) ysc(on) xlabel(0 .0005 .001 .0015 .002 .0025 .003) xtitle("Probability of Dispute Initiation") ylabel(0 20 40 60) ytitle(Frequency) scheme(s1mono) 
{res}{txt}
{com}.         capture drop b1-b12 nodef def def_change
{txt}
{com}.         
. * 4
. 
.         probit dispute ptargdef [iweight=cem_weights_4]

{res}{txt}Iteration 0:{space 3}log likelihood = {res:-8770.9635}  
Iteration 1:{space 3}log likelihood = {res:  -8674.91}  
Iteration 2:{space 3}log likelihood = {res:-8672.6408}  
Iteration 3:{space 3}log likelihood = {res:-8672.6385}  
Iteration 4:{space 3}log likelihood = {res:-8672.6385}  
{res}
{txt}Probit regression{col 49}Number of obs{col 67}= {res}   546,550
{txt}{col 49}LR chi2({res}1{txt}){col 67}= {res}    196.65
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log likelihood = {res}-8672.6385{txt}{col 49}Pseudo R2{col 67}= {res}    0.0112

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}     dispute{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}ptargdef {c |}{col 14}{res}{space 2}-.2699214{col 26}{space 2}  .020019{col 37}{space 1}  -13.48{col 46}{space 3}0.000{col 54}{space 4}-.3091579{col 67}{space 3}-.2306849
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-2.740257{col 26}{space 2} .0108709{col 37}{space 1} -252.07{col 46}{space 3}0.000{col 54}{space 4}-2.761563{col 67}{space 3} -2.71895
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}.         estsimp probit dispute ptargdef pchalally pchaloff pchalneu ln_distance capprop jdem s_un_glo peaceyrs peaceyrs2 peaceyrs3 [iweight=cem_weights_4]

{txt}Iteration 0:   log likelihood = {res}-8770.9635
{txt}Iteration 1:   log likelihood = {res}-7917.5967
{txt}Iteration 2:   log likelihood = {res}-7718.1513
{txt}Iteration 3:   log likelihood = {res}-7710.5722
{txt}Iteration 4:   log likelihood = {res}-7710.5423
{txt}Iteration 5:   log likelihood = {res}-7710.5423

{txt}Probit regression                                 Number of obs   = {res}    546550
                                                  {txt}LR chi2({res}11{txt})     = {res}   2120.84
                                                  {txt}Prob > chi2     = {res}    0.0000
{txt}Log likelihood = {res}-7710.5423                       {txt}Pseudo R2       = {res}    0.1209

{txt}{hline 13}{c TT}{hline 64}
     dispute {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
    ptargdef {c |}  {res}-.2992295   .0217482   -13.76   0.000    -.3418551   -.2566039
   {txt}pchalally {c |}  {res}-.4329878   .0231263   -18.72   0.000    -.4783144   -.3876611
    {txt}pchaloff {c |}  {res} .3353802   .0736823     4.55   0.000     .1909655    .4797949
    {txt}pchalneu {c |}  {res} .1975862   .0639036     3.09   0.002     .0723374    .3228349
 {txt}ln_distance {c |}  {res}-.4480311   .0128762   -34.80   0.000    -.4732681   -.4227941
     {txt}capprop {c |}  {res} .2478262   .0283106     8.75   0.000     .1923385    .3033139
        {txt}jdem {c |}  {res}-.2326032   .0667047    -3.49   0.000    -.3633421   -.1018643
    {txt}s_un_glo {c |}  {res} -1.44285   .0777947   -18.55   0.000    -1.595325   -1.290376
    {txt}peaceyrs {c |}  {res}-.0333091   .0023775   -14.01   0.000    -.0379688   -.0286493
   {txt}peaceyrs2 {c |}  {res} .0006521   .0000673     9.68   0.000     .0005202    .0007841
   {txt}peaceyrs3 {c |}  {res}-3.31e-06   4.67e-07    -7.09   0.000    -4.23e-06   -2.40e-06
       {txt}_cons {c |}  {res} 2.351282   .1456122    16.15   0.000     2.065887    2.636677
{txt}{hline 13}{c BT}{hline 64}
Note: 3 failures and 0 successes completely determined.

{res}Simulating main parameters.  Please wait....
% of simulations completed: 8% 16% 25% 33% 41% 50% 58% 66% 75% 83% 91% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12
{txt}
{com}. 
.         setx mean if ptargdef==0
{txt}
{com}.         setx ptargdef 0
{txt}
{com}. 
.         simqi, prval(1) genpr (nodef)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0016142     .0000805     .0014606     .001777

Simqi generated the following new variable(s): nodef
{txt}
{com}.         
.         _pctile nodef, p(2.5,97.5)
{txt}
{com}.         replace nodef =. if nodef < r(r1) | nodef > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize nodef

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}nodef {c |}{res}        950    .0016141    .0000704   .0014612   .0017768
{txt}
{com}. 
.         setx mean if ptargdef==1
{txt}
{com}.         setx ptargdef 1
{txt}
{com}. 
.         simqi, prval(1) genpr (def)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0005892     .0000432     .0005083    .0006724

Simqi generated the following new variable(s): def
{txt}
{com}.         
.         _pctile def, p(2.5,97.5)
{txt}
{com}.         replace def =. if def < r(r1) | def > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize def

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 9}def {c |}{res}        950    .0005888    .0000375   .0005084   .0006718
{txt}
{com}. 
.         gen def_change = ((def-nodef)/nodef)*100
{txt}(1,077,091 missing values generated)

{com}.         sum def_change

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 2}def_change {c |}{res}        901   -63.48838    2.457584  -70.46239  -56.08044
{txt}
{com}.         
.         twoway hist def, frequency color(gray 8) ||hist nodef, frequency color(black) legend(off) ysc(on) xlabel(0 .0005 .001 .0015 .002 .0025 .003) xtitle("Probability of Dispute Initiation") ylabel(0 20 40 60) ytitle(Frequency) scheme(s1mono) 
{res}{txt}
{com}.         capture drop b1-b12 nodef def def_change
{txt}
{com}.         
. * We estimate the JL model using the KVP control variables (Table A7 & Figure A3)
. 
.         estsimp probit dispute ptargdef pchaloff pchalneu ptargoff ptargneu contig c_majpow t_majpow jdem rivalyears midcount  

{txt}Iteration 0:   log likelihood = {res}-16773.772
{txt}Iteration 1:   log likelihood = {res}-15440.708
{txt}Iteration 2:   log likelihood = {res}-12197.596
{txt}Iteration 3:   log likelihood = {res}-12040.572
{txt}Iteration 4:   log likelihood = {res}-12038.837
{txt}Iteration 5:   log likelihood = {res}-12038.836

{txt}Probit regression                                 Number of obs   = {res}   1077992
                                                  {txt}LR chi2({res}11{txt})     = {res}   9469.87
                                                  {txt}Prob > chi2     = {res}    0.0000
{txt}Log likelihood = {res}-12038.836                       {txt}Pseudo R2       = {res}    0.2823

{txt}{hline 13}{c TT}{hline 64}
     dispute {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
    ptargdef {c |}  {res}-.0409999   .0164652    -2.49   0.013     -.073271   -.0087287
    {txt}pchaloff {c |}  {res} .0794489   .0326763     2.43   0.015     .0154045    .1434934
    {txt}pchalneu {c |}  {res} .0202673   .0295821     0.69   0.493    -.0377125    .0782471
    {txt}ptargoff {c |}  {res} .1016525   .0360848     2.82   0.005     .0309277    .1723774
    {txt}ptargneu {c |}  {res}-.0642607   .0347403    -1.85   0.064    -.1323505    .0038291
      {txt}contig {c |}  {res} .8382703    .019677    42.60   0.000     .7997041    .8768365
    {txt}c_majpow {c |}  {res} .5450971   .0196981    27.67   0.000     .5064895    .5837047
    {txt}t_majpow {c |}  {res} .2828838   .0230568    12.27   0.000     .2376932    .3280744
        {txt}jdem {c |}  {res}-.0375124   .0270789    -1.39   0.166    -.0905861    .0155613
  {txt}rivalyears {c |}  {res}  .112049   .0051831    21.62   0.000     .1018903    .1222077
    {txt}midcount {c |}  {res} .5446933    .012522    43.50   0.000     .5201507     .569236
       {txt}_cons {c |}  {res}-3.218296   .0143803  -223.80   0.000    -3.246481   -3.190111
{txt}{hline 13}{c BT}{hline 64}

{res}Simulating main parameters.  Please wait....
% of simulations completed: 8% 16% 25% 33% 41% 50% 58% 66% 75% 83% 91% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12
{txt}
{com}. 
.         setx mean if ptargdef==0
{txt}
{com}.         setx ptargdef 0
{txt}
{com}. 
.         simqi, prval(1) genpr (nodef)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res} .0009166     .0000408     .0008355    .0010011

Simqi generated the following new variable(s): nodef
{txt}
{com}.         
.         _pctile nodef, p(2.5,97.5)
{txt}
{com}.         replace nodef =. if nodef < r(r1) | nodef > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize nodef

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 7}nodef {c |}{res}        950    .0009165    .0000355   .0008357    .001001
{txt}
{com}. 
.         setx mean if ptargdef==1
{txt}
{com}.         setx ptargdef 1
{txt}
{com}. 
.         simqi, prval(1) genpr (def)

{txt}      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
             Pr(dispute=1) |  {res}  .000738     .0000327     .0006755    .0008029

Simqi generated the following new variable(s): def
{txt}
{com}.         
.         _pctile def, p(2.5,97.5)
{txt}
{com}.         replace def =. if def < r(r1) | def > r(r2)
{txt}(50 real changes made, 50 to missing)

{com}.         summarize def

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 9}def {c |}{res}        950    .0007379    .0000286   .0006755   .0008028
{txt}
{com}. 
.         gen def_change = ((def-nodef)/nodef)*100
{txt}(1,077,089 missing values generated)

{com}.         sum def_change

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 2}def_change {c |}{res}        903   -19.38411    4.062536  -30.31806  -7.020514
{txt}
{com}.         
.         twoway hist def, frequency color(gray 8) ||hist nodef, frequency color(black) legend(order(1 "Target with Defense Pact" 2 "Target without Defense Pact")) ysc(on) xlabel(0 .0005 .001) xtitle("Probability of Dispute Initiation") ylabel(0 20 40 60) ytitle(Frequency) scheme(s1mono) 
{res}{txt}
{com}.         capture drop b1-b12 nodef def def_change
{txt}
{com}.         
. exit

{txt}end of do-file


{com}. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}C:\Users\jcjo228\Dropbox\DeterResponse\LeedsJohnsonJOPrep\LeedsJohnsonJOPrep.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res}16 Aug 2015, 15:54:08
{txt}{.-}
{smcl}
{txt}{sf}{ul off}